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Abstract 

The Fermi Large Area Telescope (LAT) collaboration recently released the updated results of the 
measurement of the cosmic ray electron (CRE) spectrum and published its first constraints on 
the CRE anisotropy. With respect to the previous Fermi-LAT results, the CRE spectrum mea- 
surement was extended down from 20 to 7 GeV, thus providing a better lever arm to discriminate 
theoretical models. Here we show that the new data strengthen the evidence for the presence 
of two distinct electron and positron spectral components. Furthermore, we show that under 
such hypothesis most relevant CRE and positron data sets are remarkably well reproduced. Con- 
sistent fits of cosmic -ray nuclei and antiproton data, which are crucial to validate the adopted 
propagation setup(s) and to fix the solar modulation potential, are obtained for the Kraichnan 
and plain-diffusion propagation setups, while the Kolmogorov one is disfavored. We then con- 
firm that nearby pulsars are viable source candidates of the required e ± extra-component. In 
that case, we show that the predicted CRE anisotropy is compatible with Fermi-LAT constraints 
and that a positive detection should be at hand of that observatory. Models assuming that only 
nearby supernova remnants contribute to the high energy tail of the observed CRE spectrum are 
in contrast with anisotropy limits. 



1. Introduction 

Last year, the Fermi-LAT Collaboration published the electron + positron spectrum in the 
energy range between 20 GeV and 1 Te V, measured during the first six months of the Fermi mis- 
sion 111]. That result came in the middle of a hot debate which arose as ATIC |2] and PAMELA 
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1 3] collaborations reported some anomalies in the Cosmic Ray (CR) leptonic sector. In particu- 
lar, ATIC observed a pronounced bump in the electron + positron spectrum at around 600 GeV, 
while PAMELA found the positron fraction e + /(e + + e~) to increase with energy above 10 GeV. 
Both features are hardly compatible with the standard scenario in which CR electrons (e~) are 
accelerated in the surrounding of Galactic supernova remnants (SNRs) and positrons are pre- 
dominantly of secondary origin (i.e. they are produced only by spallation of CR nuclei onto the 
interstellar medium gas). Below 100 GeV the Cosmic Ray Electron (CRE) spectrum observed 
by the Fermi-LAT agrees with the one measured by ATIC, but at higher energies it does not dis- 
play the feature seen by that experiment, being compatible with a single power law with index 
3.045 + 0.008. The absence of the ATIC feature was confirmed by the H.E.S.S. atmospheric 
Cherenkov telescope [4]. Below 1 TeV the spectrum measured by H.E.S.S. is in agreement with 
Fermi-LAT's while above that energy (a range not probed by Fermi-LAT) a significant spectral 
steepening was reported by the former experiment. 

Soon after that measurement was released, several members of the Fermi-LAT Collaboration 
showed - using the GALPR0P package 0] - that an interpretation of the observed spectra is 
possible within a conventional model in which SNRs are the only primary sources of CRE [7]. 
Below 20 GeV, however, they found the predictions of that model to be in tension with pre-Fermi 
measurements, in particular HEAT [8] and AMS-01 [9] data. Furthermore, the positron fraction 
measured by PAMELA was not reproduced within that framework neither above (where its ob- 
served slope is opposite to the theoretically expected one) nor below ~ 10 GeV. An alternative 
scenario was therefore proposed in the same paper J7[], which invoked the presence of an extra 
component of electrons and positrons. In order to reproduce the observed CRE spectrum and 
positron fraction, the injection spectrum of such extra component has to be harder than the one 
expected to be accelerated by SNRs and must display a cutoff at around 1 TeV. In that case it 
was shown that both Fermi-LAT and PAMELA measurements above 20 GeV could be explained 
consistently and the disagreement of the model with pre-Fermi data below 20 GeV was reduced 
(the remaining discrepancy was guessed to be a consequence of solar modulation). 

Following the suggestion of some pioneering works iTToL 11] (see also the more recent [ 12]), 
in J3] it was argued that such extra-component can naturally be produced by middle-aged pulsars 
in the vicinity of the Solar System. An alternative scenario based on the annihilation of dark 
matter particles with mass around 1 TeV/c 2 was also discussed in the same paper and shown 
to be compatible with Fermi-LAT, H.E.S.S. and PAMELA data (see also II 1311 and a number of 
following papers, e.g. [l5lo . Several counter-arguments, however, have been risen which 
make this interpretation disfavored if not ruled out (see e.g. [16]). Other possible origins of 
the required electron and positron extra-component have also been proposed, e.g., enhanced 
secondary production in standard SNRs [17] or an inhomogeneous distribution of source in the 
Solar System neighborhood lHHl . 

Recently, the Fermi-LAT Collaboration released a new measurement of the CRE spectrum 
based on one year data. The observed spectrum extends down to 7 GeV lfl9ll and is confirmed to 
be compatible with a single power-law with a spectral index 3.08 + 0.05, only slightly softer than 
the one published in [1]. Hints of a deviation from a pure power-law behavior between 20-100 
GeV and at a few hundreds GeV, which were found in the six month data, are still present in 
the updated spectrum. The most relevant new piece of information comes from the data in the 7 
- 20 GeV energy range. Interestingly, the updated Fermi-LAT spectrum agrees with HEAT [8] 
and AMS-01 [9] data. As we will show these findings rise serious troubles to single component 
models and call for some revisions also for the two-component model discussed in 0]. 

Such problems become even more challenging if one considers complementary CR data sets, 
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in particular the positron-to-electron fraction and the antiproton spectrum measured by several 
other experiments. Interestingly, even at low energy several discrepancies can be found between 
the data and the predictions of numerical CR diffusion models, as GALPROP or DRAGON Q 1 2(J, if 
a conventional Kolmogorov propagation setup is adopted. Indeed, the positron fraction measured 
by PAMELA is significantly lower than the numerical estimates obtained within this setup below 
10 GeV (see e.g. |13j,|7[]). A long-standing problem is represented by the antiproton deficit which 
was found by the GALPROP team with respect to several experimental data J2lll . Recently it 



became more serious with the release of PAMELA antiproton high statistics data 11221 12311 (see 



urn- 

A possible solution of some of those anomalies has been recently proposed in terms of 
charge-dependent solar modulation [26]. To explain Fermi-LAT data, however, this scenario 
requires to assume quite extreme conditions during the latest solar activity phase. Furthermore, 
its predictions disagree with the preliminary measurements of the absolute electron (e~) spectrum 
performed by PAMELA below ~ 7 GeV 12711 . It is important to mention here that Fermi-LAT 
and PAMELA have been taking data during the same solar cycle, which allows to reduce the 
uncertainties related to solar modulation. 

Here we aim at obtaining a consistent description of all the observations mentioned above 
following a different approach: we work in a simple force field charge independent modulation 
framework but we consider other propagation setups besides the Kolmogorov one adopted in 
flEii- In particular, we adopt a model based on Rraichnan-like turbulence with a lower level 



of reacceleration with respect to the Kolmogorov one: this model was presented in [24] and 
obtained as a result of a combined maximum likelihood analysis based on B/C data (including 
the recent CREAM data-set B28IO and p/p updated measurements. We will show that a consistent 
description of most CR data sets is indeed possible for models of this kind. 

This work is organized in the following way: in Sec. [2] we will identify three reference 
propagation setups compatible with cosmic -ray nuclei observations; then in Sec. [3] we will adopt 
those setups trying to build single-component CRE models and showing why they do not provide 
a satisfactory description of the available data; in Sec.|4]we will describe how the introduction of 
an e ± extra-component allow a consistent interpretation of those data; finally, in Sec. [5] we will 
show how pulsars can rather naturally power such extra-component, and check if the proposed 
scenario is compatible with upper limits on the CRE anisotropy recently released by the Fermi- 
LAT collaboration 12911 . Sections 6 and 7 are devoted to a discussion of our findings and to 
concluding remarks. 

2. The propagation setups 

The CR propagation equation from a continuous distribution of sources can be written in the 
general form: 

= Q\p, r,Z) + J] c yS« g as(r, z)o- ft N j - c^cr^E^N 1 . ( 1 ) 



'Code available at http://www.desy.de/~maccione/DRAGON/ 
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Here N'(p, r, z) is the number density of the z'-th atomic species; p is its momentum; p 1 the particle 
velocity in units of the speed of light c; cr,„ is the total inelastic cross section onto the inter 
stellar medium (ISM) gas, whose density is n gas ; cr^ is the production cross-section of a nuclear 
species j by the fragmentation of the z'-th one; D is the spatial diffusion coefficient; v c is the 
convection velocity. The last term on the l.h.s. of Eq. (Q]) describes diffusive reacceleration of 
CRs in the turbulent galactic magnetic field. In agreement with the quasi-linear theory we assume 
the diffusion coefficient in momentum space D pp to be related to the spatial diffusion coefficient 
by the relationship (see e.g. IBOIO 

D - A 1 (2) 

pp 36(4 - £ 2 )(4 - 6) D 

where Va is the Alfven velocity and 5 the index of the power-law dependence of the spatial diffu- 
sion coefficient on the particle rigidity (see Eqj4|. Here we assume that diffusive reacceleration 
takes place in the entire diffusive halo. 

The distribution of Galactic CR sources is poorly known. The large scale source distribution 
which we adopt here is approximated to be cylindrically symmetric, so that the source term for 
every nuclear species i takes the form: 

Qi(E k , r, z) = fs (r, z) goj " , (3) 

For CR nuclei, we assume that fs (r, z) traces the distribution of supernova remnants and impose 



the normalization condition fs (r , z©) = 1 . We use the same distribution as in 12411 for all models 
considered in this paper. This is slightly different from that used in GALPROP which, however, 
does not have significant effects on the CRE and positron spectra. We also assume that the 
injection spectral index is independent of the primary nucleus and is the same as that of protons 
y, — y p . Since SNRs are expected to accelerate CR nuclei up to at least 10 15 eV, which is well 
above the energies we consider here, we do not account for a possible high energy suppression 
of the injection spectrum. 

We neglect CR convection as we tested that under reasonable conditions its effects are neg- 
ligible at the energies under consideration here and that it is not necessary to reproduce experi- 



mental data. Nuclei spallation is treated as described in 161 13 1L 13211 . 

Here we assume the diffusion coefficient D to be spatially uniform and that it only depends 
on the particle rigidity p and the particle speed f3 according to the following relation 

D(p) = D /?>(—) , (4) 



Po 

with 77 controlling essentially the low energy behavior of D. While one would expect 77 = 1 as 
the most natural dependence of diffusion on the particle speed, it should be taken into account 
that diffusion may actually be inhibited at low energies due to the back-reaction of CR on the 
magneto-hydrodynamic waves. In a dedicated analysis [33] of that effect, a low-energy increase 
of D was found. While such a behavior cannot be represented as a simple function of /3 and p, 
an effective value of 77 can, nevertheless, be found which allows to fit low energy data. Clearly, 
the precise value depends on the details of the model. For example in [33] 77 — -3 was found to 
fit the data, while the authors of ll34Tl found 77 = -1.3, in both cases for models with 5 = 0.5 (but 
rather different values of other parameters). In i24ll we found that 77 = -0.4 allows a rather good 
fit of low energy nuclear data for models with low reacceleration and 6 0.5. 
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Table 1 : Propagation and CR injection parameters for the three reference models considered in this paper. Dq is the 
diffusion coefficient normalization at 3 GV; 6 the index of the power-law dependence of D on energy; z.h the half-height 
of the Galactic CR confinement halo; y p the CR nuclei injection index; va is the Alfven velocity; r\ is the exponent of the 
power-law dependence of D on the particle velocity /}. The source spectral break rigidity for the KOL model is 4 GV. 
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KRA 


0.5 


3.0 x 10 28 
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2.25 


v A = 15 


-0.4 


PD 


0.60 


2.4 x 10 28 


4 


2.15 


V A =0 


-0.4 



Although the description of propagation offered by Eq. (Q} is already extremely simplified 
and in fact just an effective one, present CR data are still not accurate enough to fix the values of 
the main parameters controlling propagation: Do and 6 from Eq. va from Eq. (0, the height 
of the Galactic diffusion region Zh, and the injection index y p appearing in Eq. (f3j). Moreover, 
when considering data below a few GeV/n also the parameter rj (see Eq. (O) and the modulation 
due to solar activity play a significant role and must be taken into account. For these reasons, here 
we consider three reference propagation models defined by the parameters reported in Tab.Q] In 
making this choice, we fixed zi, = 4 kpc.0 Models are defined as follows: the PD model is a plain 
diffusion one, in which we tried to reproduce CR data with only diffusion and spallation losses, 
hence setting va = 0; the KOL model is built assuming that the power spectrum of the galactic 
turbulent magnetic field is of Kolmogorov type, which fixes 5 — 1/3, while the KRA model 
assumes a Kraichnan spectrum for the galactic turbulent magnetic field, hence setting 6 - 1/2. 
With these choices being made for the physical models, the remaining free parameters were tuned 
to reproduce the boron to carbon ratio (see Fig.Q]) in the energy range 0. 1 < E < 10 3 GeV/n. We 
take experimental data from HEAO-3 lf35ll and CRN Ir36ll satellite-based experiments and from 
CREAM 111 and ACE/CRIS Q. HEAO-3 B/C data are nicely confirmed by a recent analysis 
of AMS-01 data IH. 

We treat the solar modulation adopting the usual force-free field approach [40] and tune 
the modulation potential <1> so as to reproduce the proton spectrum measured by PAMELA i39ll 
during the same solar cycle period during which Fermi-LAT performed the CRE spectrum mea- 
surement (see Fig.|2|. 

It is clear from Tab. Q] and from Fig. Q] that the steeper the energy dependence of D (i.e. the 
larger the value of 6) the lower is the amount of reacceleration required to reasonably match 
the B/C data. Indeed, if va =0 we need 6 = 0.6, while if va = 30 km/s then 6 = 0.33 is 
enough to reproduce high energy B/C data. Although the interpretation of low energy data is 
not straightforward, because of several poorly known effects (e.g., convection or dissipation of 
Alfven modes) which are expected to play an increasing role with decreasing energy, we found 
that a proper tuning of the parameter n is enough to provide an effective description of CR data 
even well below the GeV (see also [24]). 

Although the B/C ratio is the most used quantity for the purpose of fixing CR propagation 
parameters, also other quantities, in particular the p/p ratio or the antiproton absolute spec- 
trum provide complementary constraints. Indeed, ratios of nuclear species like B/C are weakly 
insensitive to the injection spectra of primaries, since the secondaries have almost the same en- 
ergy/nucleon as the primaries. Hence we have still the freedom to adjust the injection spectra to 



2 We verified that, as long as 1 kpc < zi, < 10 kpc, changing zt, almost amounts to a redefinition of the value of Dq. 
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Figure I: B/C ratios, as computed with our three reference models given in Tab.fT] are compared with experimental data. 
Red lines: Kolmogorov model (KOL); blue lines: Kraichnan model (KRA); green lines: plain diffusion (PD). Solid lines 
are modulated with (D = 450 MV, which is appropriate for HEAO-3 and CREAM data, while triple-dotted lines have 
been computed for il> = 300 MV which is appropriate for ACE data. 



reproduce the proton and other primary spectra. Antiprotons, on the other hand, are produced in 
a spectrum by the spallation of the primary protons and Helium nuclei, which depends signifi- 
cantly on the energy of the primary particle. Therefore, antiprotons are sensitive to the injection 
spectrum of primaries, and break the degeneracy. 

Constraints coming from antiproton data became quite stringent with the recent release of 
PAMELA data |22[ 12311 and cannot be ignored. PAMELA proton data are also crucial to set, 
for each model, the solar modulation potential. It is interesting that all models considered here 
match those data for the same value of the modulation potential (see Fig.[2]i0. 

Once the modulation potential has been fixed in that way, model predictions can be compared 
with experimental antiproton data. It is clear from Fig. [3] that while the propagation models 
considered in the above are almost degenerate against the B/C data some of them, in particular 
the KOL model, are disfavored by the antiproton data. This was also established on the basis of 
a recent combined statistical analysis of CR nuclei and antiproton data B2411 . 

Overall, this seems to favor models with low reacceleration and a value of the diffusion 
coefficient slope close to 0.5 or higher, which is also in agreement with independent findings 



1 2.511 - In the next section we will see as CR electron and positron data also lead to a similar 



conclusion. 



3. Modeling the CRE spectrum with a single Galactic component 

We start our analysis by trying to interpret CR electron and positron data with single-component 
models. 



3 It should be noted that differently from the KRA and PD models, the KOL model requires a pronounced break in the 
proton injection spectrum in order to reproduce the data (see Tab. 1 and Fig.f2j- 
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Figure 2: The proton spectra calculated for our three reference models. Color notation is the same as in the previous 
figure. The solar modulation potential, which was tuned to reproduce experimental during solar activity minimum, is 
d> = 550 MV for all models. Dashed lines are the corresponding local interstellar (LIS) spectra. 
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Figure 3: The theoretical antiproton spectra computed for our three reference models are compared with experimental 
data. The color notation is the same as in the previous figures. All solid lines are modulated according to the force field 
scheme assuming the charge symmetric potential <I> = 550 MV as required to reproduce the measured proton spectrum 
at low energies (see Fig. [2}- Dashed lines are the corresponding LIS spectra. 
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We evaluated these model with the DRAGON numerical package lf24Tl . We also verified that 
our results are reproduced by the GALPROP package [5J |6J] under the same physical conditions. 
We consider the three diffusion setups, KOL, KRA and PD discussed in the previous section and 
for each of them we tune the free parameters involved in the calculation of the CRE spectrum 
against Fermi-LAT data in the energy range 20-100 GeV, where both statistics and systematics 
errors are smallest. The choice to restrict our model fitting to data in this energy range is also 
strongly motivated by theoretical arguments. Indeed, above 20 GeV the uncertainties on the CRE 
spectrum due to solar modulation are negligible; below 100 GeV the CRE spectrum is hardly 
affected by the uncertainties due to the stochastic spatial and temporal distribution of CRE local 
sources which instead become relevant in the TeV energy range I4ll l7ll. 

We find that above 20 GeV the injection spectral indexes y - 2.45, 2.37, 2.32, respectively 
for the KOL, KRA and PD models, are adequate to reproduce Fermi-LAT data (see Fig. |4(a)| >. 

Problems appear, however, when considering lower energy data. A discrepancy between 
low energy pre-Fermi data and the prediction of single component models was already noticed 
in J3] and tentatively ascribed to solar modulation or to systematics experimental errors. That 
interpretation, however, does not hold against the new Fermi-LAT data between 7 and 20 GeV 
which agree with AMS-01 and HEAT data and are clearly incompatible with the results of the 
above models. This is the case both if single power-law source spectra are adopted and if a 
spectral break is introduced, as shown in Fig.s |4(a)] and |4(b)| respectively. 

It is worth noticing that reacceleration models need spectral breaks to correct for the anoma- 
lous behavior which would otherwise arise in the propagated LIS spectrum. The reason of such 
a behavior can be traced back to the combined effect of reacceleration and energy losses. Reac- 
celeration shifts electrons from the low-energy to the high-energy region of the spectrum, while 
energy losses take electrons from high to low energy. The two effects have comparable strength 
in the GeV region and give rise to pronounced bumps in the unmodulated spectra shown in 
Fig. |4(a)| if a single power-law is assumed. Clearly, this feature is more evident in models with 
strong reacceleration (as in the KOL model), and must be treated by introducing a sharp, and 
hardly justifiable, break in the injection. On the other hand, the KRA model, in which only mod- 
erate reacceleration is present, requires a smoother break. No break at all is required for the PD 
model. 

As it is evident from Fig.s |4(a)] and |4(b)} a very small modulation potential (smaller than ~ 
100 MV) would allow to reproduce low energy Fermi-LAT data. Two major problems, however, 
arise against this option: a) such a low modulation potential is not compatible with proton and 
nuclei data; b) it does not fit AMS-01 data which were taken during a solar minimum phase like 
Fermi-LAT ones were. 

In principle, both problems could be circumvented in a charge asymmetric modulation sce- 
nario, in which solar modulation is expected to act differently on particles of opposite electric 
charge during periods of opposite polarity of the Solar magnetic field. This is the case for AMS- 
01 and Fermi data taking periods. Interestingly, the antiproton and positron fraction measured 
by PAMELA and other experiments can be consistently reproduced in that framework [26]. This 
possibility, however, seems to be at odds with recent (still preliminary, see M27I0 measurements 
of the electron e~ spectrum by PAMELA, which agree both with Fermi-LAT, between 7 and 100 
GeV, and with AMS-01 below 10 GeV. Indeed, even adopting a vanishing potential for the e~ 
(which is a quite extreme assumption), such models fail to match PAMELA e~ data down to a 
few GeV. Therefore, since PAMELA and Fermi-LAT data were taken during the same solar cycle 
phase, it is evident that the discrepancy between the prediction of single component models and 
low energy data cannot be ascribed to charge asymmetric solar modulation. 
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Figure 4: The e + e + spectra computed for single component models with unique injection slopes y = 2.45, 2.37, 2.32 
(panel a) and broken power-law y = 1.6/2.4; 2.0/2.37; 2.32/2.32 below/above 2 GeV for the KOL (red), KRA (blue) 
and PD (green) diffusion setups respectively (panel b). Fermi-LAT data are represents by the red points. See Fig. 5a,c,e 
for the other experiment's symbols. Modulated lines (solid curves) have been computed in the force field framework 
with <D = 550 MV; dashed lines are the corresponding local interstellar (LIS) spectra. 
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A reasonable fit of the CRE spectrum measured by Fermi-LAT is possible only at the price of 
normalizing the models to data at 10 GeV (rather than 100 GeV) and adopting injection spectra 
slightly steeper than those reported in the above (2.50 rather than 2.45 in the KRA model, 2.43 
rather than 2.37 for the KOL model). However, this is not the most natural choice for the reasons 
we explained above. Moreover, the observed spectral slope between 7 and 100 GeV is never 
reproduced in those cases, as it is clear from Fig. |5(a)[|5(c)| and |5(e)] >. 

For the same models, we also computed the corresponding positron fraction e + /(e~ + e + ) (see 
Fig. [5(5)1 |5(d)l and[5^)). Clearly, below 10 GeV the e + /(e~ + e + ) measured by PAMELA can be 
reproduced by the KRA and PD models while the fit is unsatisfactory for the KOL model. Again, 
low reacceleration models seem to provide a better description of low energy experimental data. 
None of the single component realizations, however, can reproduce PAMELA data above 10 
GeV. 

We conclude this section by summarizing the main drawbacks of single component CRE 
models: 

• they are unable to reproduce all the features revealed by Fermi-LAT in the CRE spectrum, 
in particular the flattening observed at around 20 GeV (which was also recently found 
by PAMELA 112711 ) and the softening at ~ 500 GeV. If, as they should, such models are 
normalized against data in an energy range where systematical and theoretical uncertainties 
are the smallest, they clearly fail to match CRE Fermi-LAT and PAMELA data below 20 
GeV. Most importantly, the CRE spectral slope between 7 and 100 cannot be reproduced. 

• They suffer from relevant problems also at very high energy, because they do not repro- 
duce the softening of the CRE spectrum measured by HE.S.S. above 1 TeV This requires 
either to introduce a cutoff in the CRE source spectrum, which, however, does not account 
for the observed spectral hardening around 100 GeV, or to introduce an additional CRE 
component of local origin, with a spectrum peaked around the TeV, as first suggested in 



111]. 



As explained in many papers (see e.g. 14211 ). they cannot reproduce the rising positron-to- 
electron ratio measured by PAMELA at high energy. Therefore, if PAMELA observations 
are correct, an additional positron component, besides that produced by CR spallation, has 
to be invoked. 



4. Two components models 

In the following we try to reproduce Fermi-LAT, HE.S.S. and PAMELA CRE and positron 
data in the framework of a two component scenario. The approach we follow in this section is a 
rather straightforward generalization of what we did in the previous one: we assume the presence 
in the Galaxy of two CRE source populations rather than one. The first component accelerates 
only electrons (the standard component) with a spectral slope y e - > 2 (as expected from Fermi 
acceleration in SNRs) up to some cutoff E e cuV while the second one accelerates electrons and 
positrons with a common spectrum (extra-component) with a form 

Qe<E) = Q ^-\ * e- E/E ^ , (5) 

and y e ± < 2 as required to account for the positron fraction rise observed by PAMELA. Here 
we also assume that both source classes have the same continuous spatial distribution as that 
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Figure 5: The electron (g~ + e + ) spectrum is shown for the KOL (panel a), KRA (panel c) and PD (panel e) diffusion 
setups as specified in Tab.^ The electron source spectral indexes are 1.60/2.50 below/above 4 GeV for the KOL model 
(panel a) and 2.00/2.43 and 2.0/2.400 below/above 2 GeV for the KRA (panel c) and PD models (panel e) respectively. 
Purple and orange lines represent secondary electron and positron spectra respectively. The corresponding positron 
fraction (e + l(e~ + e + )) curve computed under the same conditions are shown in the panels (b), (d) and (f) respectively. 
Solid lines are modulated with a potential <I> = 550 MV. Dashed lines are the corresponding LIS spectra. 
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considered in the previous section. This assumption, however, is not critical, and we checked 
that a quite wide sample of extra component source distributions (including the one expected 
from the annihilation of dark matter with a NFW profile) yields quite similar results after a 
small tuning of the injection parameters. The assumptions of considering the source distribution 
continuous and cylindrically symmetric up to ~ 1 TeV are more crucial in this respect and will 
be reconsidered in the next Sections. 

In order to tune the normalization of the two spectral components we follow a multi-messenger 
approach. Our first step is to tune the standard component to consistently reproduce both the 
e~ +e + spectrum measured by Fermi-LAT and the e + /(e~ + e + ) measured by PAMELA below 20 
GeV, where the effect of the extra component is supposed to be negligible. Remarkably, this is 
indeed possible if we use propagation setups with low reacceleration, namely either the KRA or 



the PD setups. Noticeably, these setups (especially the KRA one [24]) also provide the best com- 
bined fits of the B/C data, proton and the antiproton spectra measured by PAMELA (see Sec. [2]). 
Fermi-LAT and PAMELA electron, positron fraction and antiproton data are all reproduced with 
the same modulation potential <S> = 550 MV, in the simple force field framework. The required 
source spectral slopes for the electron standard component is y e - = 2.00/2.65, 2.00/2.60 be- 
low/above 4 GeV for the KRA and PD setups respectively. We also introduce an exponential 
cutoff in the source spectrum of this component at 3 TeV. The choice of higher cutoff would not 
affect significantly or final results. 



Clearly, in the absence of the extra e ± component, high energy CRE (Fig. 6(a)[ ) and positron 



fraction (Fig. |6(b j| > data would completely be missed (see dotted line in Fig. |6(a)| >. Remarkably, 
the indication of the presence of an extra CRE component could be found even considering only 
Fermi-LAT and the low energy PAMELA data alone. This a new and interesting result which 
was made possible by the spectacular data collected by Fermi-LAT and PAMELA and by the fact 
they operate during the same solar activity phase (hence reducing the uncertainties due to solar 
modulation). 

Our next step is then to suitably choose the extra-component as to reproduce Fermi-LAT 
and H.E.S.S. high energy CRE data. We find here that this is possible by taking y e ± = 1.5 and 
£cut = 1.0 H- 1.5 TeV both for the KRA setup (see Fig. |6(a)| i and the PD one. This is similar 
to what done in 0, l55ll but for the choice of the propagation setup which in those papers was 
assumed to obey Kolmogorov diffusion (as a consequence, low energy PAMELA data were not 
reproduced in that case). 

It is interesting that preliminary e~ spectrum measured by PAMELA 12711 is also nicely re- 
produced by our extra-component models. Above 100 GeV this spectrum is softer than the 
<?~ + e + measured by Fermi-LAT by the exact amount which is required to leave room the e + 
extra-component. It is remarkable that such a relatively simple approach allows to reproduce 
such a large number of observations. 

The fact that the high energy tail of the CRE spectrum is reproduced so well in terms of a 
continuous extra component source distribution suggest that either it is dominated by a physically 
continuous distribution of sources (as one expects for dark matter) or by a single astrophysical 
source in the nearby. In the next section we will discuss the latter case. 

5. The role of astrophysical nearby sources 

The nature of the extra-component of primary electrons and positrons that we invoked in the 
previous section is an intriguing matter of debate, and the possible scenarios include both an 
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Figure 6: The e~~ + e + total spectrum (panel a) and positron fraction (panel b) data are compared with the predictions 
of our two-component model. The dotted line represents the propagated standard component with injection slope y e - = 
2.00/2.65 above/below 4 GeV and E' ut = 3 TeV, while the dot-dashed line is the e ± component with y e ± = 1.5 and E' M = 
1.4 TeV. Both components are propagated with DRAGON adopting the KRA setup. The blue solid/dashed lines represent 
the modulated/LIS quantities. The modulated lines have been computed adopting the charge symmetric modulation 
potential O = 550 MV. 
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exotic explanation (involving annihilation or decay of Particle Dark Matter) or a purely astro- 
physical interpretation. Here we concentrate on the second possibility and show that well known 
astrophysical sources can rather naturally account for the available CRE electron and positron 
fraction observations. In Sec.|6]we will show that CRE anisotropy measurements should soon be 
able to confirm this possibility. 

Differently from the previous section, here we treat the extra component as originating from 
a discrete collection of sources rather than from a continuous distribution. This is the proper ap- 
proach to be followed for energies exceeding few hundreds GeV at which the loss length due to 
synchrotron emission and IC scattering becomes comparable to the average SNR mutual distance 
so that only few sources within few hundred parsecs are expected to dominate the CRE flux. To 
be even more realistic, we work with actually observed astrophysical CRE source candidates, 
pulsars and SNRs, as taken from radio catalogues. We then treat electron and positron propa- 
gation from those source to the Solar system by solving analytically the diffusion-loss equation 
similarly to what done in [7]. Since the analytical expressions and tools used here are the same 
followed in that work we do not report here any detail and address the interested reader to that 
paper. 

Concerning the large scale (standard) component, we model it with DRAGON as done in the 
previous section. For consistency, we treat analytical and numerical propagation under the same 
physical conditions. 

Since the extra component does not affect the low energy tail of the CRE and positron spectra, 
all we did in the previous section to fix the standard component holds also here. Therefore we 
again assume y e - = 2.65 and a KRA propagation setup. Similar results can be obtained with a 
PD setup. This is different from what done in [7] where the KOL setup was used. 

5.7. The contribution of local pulsars 

Pulsars are extreme astrophysical environments that release very large amounts of energy 
(~ 10 52 -J- 10 54 erg) during their active lifetime. Most of this energy is emitted in the early stage 
of their evolution, since they are powered by their spin-down. They are believed to be sources 
of both electrons and positrons (see e.g. 1 1 OL 14311 ): in fact, magnetospheric models predict the 
formation of curvature photons that - interacting with the strong pulsar magnetic field - form 
electron-positron pairs. Those particles are then expected to be accelerated at the termination 
shock of the pulsar wind nebula (PWN) and be released in the ISM when the PWN merges the 
ISM or the pulsar exit the parent SNR shell due to its proper motion (see l47ll for a detailed 
description of the latter scenario). In both cases, this process typically happens in a time between 
10 4 and 10 5 years. 

The possibility that electrons and positrons from nearby pulsars can dominate the high energy 
tail of the CRE spectrum and explain the rising behavior of the positron fraction was already 
proposed in 111 ill and studied in several more recent papers (see e.g. IU2U50I1 '). 

Similarly to the approach taken in [7], we model the CRE emission from pulsars as a point- 
like burst with a time delay with respect to the birth of the object. Such a delay is required to 
reproduce the smooth CRE spectrum measured by Fermi-LAT and H.E.S.S. and is motivated by 
the fact that electrons and positrons are expected to be trapped in the nebula until it merges with 
the ISM. We assume that electrons and positrons are accelerated in equal amounts by pulsars 
with a power-law injection energy up to an exponential cutoff at energy around 1 TeV. 

We considered the pulsars within 2 kpc from Earth, taken from the ATNF catalogue [44] Q. 



http://www.atnf.csiro.au/research/pulsar/psrcat/ 
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We verified that more distant pulsars give a negligible contribution. We also verified that y-ray 
pulsars which have been detected by Fermi-LAT 114511 and are not in the ATNF catalogue, which 
can also contribute to the observed e~ + e + spectrum (see [ 46]) do not affect significantly our 
results. 

For middle-aged pulsars, the rotational energy released at the time T can be approximated as 
. T 2 

E(T) Epsd — . where Epsd is the spin-down luminosity determined from the observed pulsar 
to 

timing and to is a characteristic braking time. Due to the strong energy losses in the PWN, only 
electron-positron pairs which are accelerated before the escape time r esc can contribute to the 

(T - T esc ) 2 

CRE flux. Therefore the e ± energy released in ISM is E e ±(T) rj e ± £psd — — , where r\ e ± 

TO 

is the e ± pair conversion efficiency. 

As shown in Fig. |7(a)| and |7(b)| the pulsar scenario allows a very good fit of both high-energy 
Fermi-LAT and H.E.S.S. electron+positron as well as PAMELA positron fraction data, due to the 
presence of primary positrons. The plots in Fig.s |7(a)| and |7(b)| have been computed assuming 
that pulsars spin down as a magnetic dipole (braking index n = 3) and that to = 10 4 yr, as 
commonly done in the literature. Furthermore we tuned the pulsar injection parameters to the 
values y e ± = 1.4, E cut = 2 TeV, r esc = 75 kyr and r] e ± = 35%. These values are compatible 
with multichannel observations of pulsars (see e.g. 1IT2I l47ll and ref.s therein). Significantly 
larger spin-down power, hence smaller values of rj e ±, are obtained for values of the braking index 
smaller than 3 as observed for several pulsars (note e.g. that for the Crab pulsar n = 2.5), which 
make this scenario even more palatable. 

While it is certainly unrealistic to assume that all pulsars share the same values of those pa- 
rameters, this is not critical for our results since the high energy tail of the spectrum is always 
dominated by a single object, namely the Monogem pulsar (PSR B0656+14). This is a conse- 
quence of its small distance (d 288^7 pc) [48], relatively young age (r = 1.1 x 10 5 yr) and 
spin-down luminosity E 3.8 x 10 34 ergs -1 . Monogem available rotational energy at the adopted 
e ± escape time is E(t > r esc ) 5 x 10 47 erg. Furthermore, it was shown in [7] that it is possible 
to reproduce CRE and positron data for several allowed combinations even randomly varying the 
pulsar parameters. 

Due to the time delay between their birth and the e ± release, very young pulsars, such as 
Vela, which are bright in the GeV and TeV gamma-ray sky, do not contribute, which explains the 
absence of pronounced bumps in the CRE spectrum around or above the TeV which were instead 
predicted in several other papers (see e.g. |49[ 50ll ). 



5.2. Adding the contribution from local SNRs 

Similar to the approach of J3l, in the above we treated the standard electron component as 
originated by a continuous distribution of sources. This approximation, however, is not realistic 
above few hundred GeV for the reasons we explained at the beginning of this Section. Indeed, if 
we share the common, and well motivated, wisdom that CRE are accelerated by SNR, only few 
of those objects will contribute to the CRE observed spectrum above 100 GeV. This may produce 
observable features in the total CRE spectrum. 

In order to study such effect, we treat CRE propagation from nearby SNRs similarly to what 
we just did for pulsars. Since the SNR lifetime is typically smaller than the propagation time, we 
consider the emission from a single SNR as a single burst simultaneous to the SNR birth. Hence, 
we consider all observed SNRs within 2 kpc as taken from the Green catalogue [51] and treat 
them as point-like <?~ sources with a power-law injection spectrum and an exponential cutoff. 
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Figure 7: The contribution from nearby pulsars (within 2 kpc) computed analytically for the KRA diffusion setup is 
added to a standard component similar to that shown in Fig. |6(a)| but the addition of an exponential cutoff at 3 TeV. Panel 
a): electron+positron spectrum. Panel b): positron fraction. Notice that since pulsars are assumed to emit both electrons 
and positrons, the rising positron-to-electron ratio measured by PAMELA is correctly reproduced within such a scenario. 
The pulsar e ± conversion efficiency is = 35% 
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Figure 8: The analytically computed CRE flux from both nearby (within 2 kpc) SNRs and pulsars is added to the same 
standard component used in Fig. |7(a)| Again, all components are propagated using the KRA diffusion-reacceleration 
setup. Panel a): electron+positron spectrum. Panel b): positron fraction. The assumed energy release for each SNR is 
taken as 2 X 10 47 erg . The pulsar efficiency is — 30%. Solar modulation modulation potential is <I> = 500MV. 



In Fig.s |8(a)| and [8(b"j| we respectively represent the CRE spectrum and positron fraction 
obtained for a reasonable combination of pulsar and SNR parameters, namely: spectral index 
y e sm = 2.4; cutoff energy Ef™ = 2 TeV; electron energy release per SN E SNR = 2 x 10 47 erg; 
T] e ± =s 30% (which is slightly smaller than that needed without considering nearby SNR) for all 
pulsars. We see from Fig. |8(a)| that under those conditions, the dominant source in the TeV 
region remains Monogem pulsars. 

Clearly, other combinations of parameters are possible, and the relative contributions of the 
several sources may vary. However, the requirement to reproduce the PAMELA positron frac- 
tion imposes an important independent constraint which does not permit to lower significantly 
the Monogem dominant contribution with respect to that of SNRs. Therefore the discrete con- 
tribution of nearby SNRs should not introduce pronounced features (bumpiness) in the CRE 
spectrum. Another relevant, and testable, prediction of this scenario is that, being the high en- 
ergy CRE flux dominated by Monogem pulsar, a significant CRE dipole anisotropy should be 
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present roughly directed towards that pulsar. This possibility will be discussed in details in the 
next Section. 

6. The CRE anisotropy 

Recently the Fermi-LAT Collaboration published upper limits on the anisotropy in electron + 
positron flux 1E9I1 . It is therefore worth checking if the models discussed in the previous section 
are compatible with such limits. 

The analytical expressions which we used to compute the CRE anisotropy due to a single 
pulsar is: 



where JV f PSR and Nf are the electron spectra from the pulsar and its sum to the large scale 
Galactic plus distant pulsar components; T and r esc are the pulsar birth time and the time it takes 
for the electrons to be released in the ISM respectively; E max (t) is the energy loss time due to 
synchrotron and IC losses (see |7|,|49|] for more details). 

The most important results we obtained are summarized in Fig. |9(a)| - 9(b)| 

First of all, it is important to notice that the model discussed in Sec. l5.1l (see Fig. |7(a)[ ), where 
only the emission from nearby pulsars is added to the smooth Galactic standard component, is 
compatible with the upper limits reported by the Fermi-LAT Collaboration. The middle-aged 
pulsar Monogem gives the dominant contribution to the anisotropy, and the expected value of 
this observable is very close to the upper limit, so that a positive detection is expected in the near 
future towards Monogem. 

Hence we now consider the model discussed in Sec. 15.21 (see Fig. |8(a)[ ), in which the main 
contribution to the high energy <?~ + e + flux still comes from pulsars but the contribution of 
local SNRs is also considered (hybrid model). We see from Fig. |9(b)| that also this possibility is 
not excluded by anisotropy measurements: in this plot, the reader may notice that the Monogem 
pulsar (red solid line) and the Vela SNR (black dashed line) contribute most to the total anisotropy 
(the black solid line), which is computed as the sum of each anisotropy weighted by the cosine 
of the angle of the corresponding source with respect to the direction of the maximum flux. 
However, also in this case the total expected anisotropy is very close to the measured upper limit, 
so that a future detection at level ~ 1 % at ~ 1 TeV towards the portion of the sky where Vela and 
Monogem are located (with the peak situated almost in the middle between the two sources) is 
to be expected in the next years. 

It should be noticed that Fermi-LAT anisotropy constraints already exclude some models 
which do reproduce the CRE and the positron fraction data. For illustrative purposes, in Fig.s 
1 10(a)| and |10(b)| we show the CRE spectrum and anisotropy computed for a a model in which 
the electron energy emitted by SNRs is increased from 2 x 10 47 erg to 1 x 10 48 erg and pulsar e ± 
emission efficiency is consequently reduced from 35% to 30% so to reproduce the observed CRE 
spectrum. We see from Fig. 1 1 0(b)] that the expected anisotropy exceeds the upper limits in that 
case, mostly because of the large contribution of Vela SNR, which is nearby (less than 300 pc) 
and very young (~ 10 4 yr). Although this result can hardly be translated into a sharp constraint 
on the e ± energy output from SNRs, due to the large number of free parameters we had to deal 
with, nevertheless we can safely conclude that a scenario in which SNRs provide the dominant 
contribution to CRE spectrum in the TeV region is incompatible with Fermi-LAT upper limits on 
the CRE anisotropy. 




(6) 
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Figure 9: Panel a): the integrated anisotropy, as a function of minimum energy, computed for the pulsar model plotted 
in Fig. |7(a)| is compared to the 95 % C.L. Fermi-LAT upper limits [29] ; Panel b) The same is done here for the pulsars 
+ SNRs model plotted in Fig. |8(a)| The black solid line represents the total anisotropy, originated mostly by Monogem 
and Vela SNR. 
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Figure 10: Panel a): The CRE spectrum for a modified version of the pulsars + SNRs (hybrid) model in which the energy 
output of SNRs is increased to 1 X 10 48 erg is represented. This model is still compatible with Fermi-LAT and HESS 
electron data. Panel b) The corresponding anisotropy is shown: a strong contribution from SNRs is not compatible with 
Fermi-LAT upper limits. 
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7. Discussion 



Here we discuss how general are the physical assumptions under which we reproduce such 
a large number of CR data sets and the compatibility of our results with other complementary 
observations. 

We start discussing the source spectrum of the electron standard component. According to 
the common wisdom we assumed that this component is accelerated by Galactic SNRs. The 
most clear evidence that electron acceleration takes place in SNRs comes from the synchrotron 
emission of those objects. The mean value of the spectral index from the Green's catalogue 
Ul is 0.5 + 0.15 in the GHz range, which loosely implies y = 2.0 + 0.3 in the 1 - 10 GeV 
range 15211 . This is therefore compatible with the source spectral index adopted in Seas [4] and [5] 
below few GeV. At larger energies, where we need y e - 2.6, we have less stringent observational 
constraints. While y-ray measurements favor harder spectral indexes above 10 GeV, it is not clear 
if the origin of that emission is leptonic. Recent y-ray observations of SNRs actually are better 
described in terms of a hadronic origin of that radiation (see e.g. Ido, 15711 '). On the theoretical 
side, Fermi linear acceleration in SNRs generally predicts values of the spectral index smaller 
than 2.5, which is marginally at odds with our finding. However, linear Fermi acceleration theory 
is likely to be not exact and the actual spectral shape could be different from that prediction. 
Furthermore, it should be noted that we modeled the standard component in the approximation of 
a cylindrically symmetric source distribution. Such description is adequate to reproduce nuclei 
spectra and secondary-to-primary ratios, but may be less realistic for electrons in the 0.1-1 
TeV range where the local distribution becomes relevant. A more realistic distribution which 
accounts for the spiral arm distribution of SNRs may actually result in a different requirement 
for the injection spectral slope fitting the data in that energy range. Indeed, being the Sun at the 
edge of a Galactic arm, the average distance from SNR is larger than in the smooth case. As 
a consequence, a harder injection spectrum is required in this case to compensate for the larger 
energy losses and reproduce the observed spectrum. We performed explorative runs probing the 
effects of a Galaxy spiral structure which actually confirm this expectation without spoiling the 
rest of our successful results. A more detailed study of this effect will be performed elsewhere. 

For what concerns the e ± extra-component, we mentioned as the hard spectrum (y e ± 1 .5) 
required to explain experimental data can be originated by pulsars wind nebulae (PWNe). While 
the acceleration mechanism responsible for such spectral shape is not understood yet, on purely 
observational grounds we know that PWNe indeed accelerate electrons with a power law spec- 
trum which is as flat as E~ l - E~ 1,s up to several hundred GeV [47]. The most critical issue here 
concerns the available energy that middle-age pulsars can release under the form of e ± pairs. We 
showed, however, that under reasonable conditions the rotational energy is sufficient to account 
for experimental data even if pairs are released in the ISM after several 10 4 yr. A physically 
viable scenario where e ± are accelerated in bow shock PWNe has been recently discussed in 
ll47Tl . The injection parameters adopted here are fully compatible with those proposed in that 
work. Forthcoming Fermi-LAT CRE spectrum and anisotropy measurements as well as y-ray 
observations of PWNe may help to validate this scenario in the next future. 

Since CR electrons contribute to the y-ray diffuse emission of the Galaxy by their bremsstrahlung 
and IC radiation, it is worth asking if the two component models presented above are compati- 
ble with Fermi-LAT observations also in that channel. For this purpose, we used our numerical 
package GAMMASKY and cross checked our results with GALPR0P. For this purpose we used the 
same radiation field as described in [59]. To derive a detailed model of the y-ray diffuse emission 
of the Galaxy is far beyond the aims of this paper. We only computed the y-ray spectrum at in- 
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termediate Galactic latitudes 10° < \b\ < 20°, averaged over all longitudes, and compared it with 
Fermi-LAT data presented in [58]. For each CR diffusion setup we computed the n component 
consistently. We treated the standard and the extra-component CRE components as described in 
Sec. [4] The contribution of local CRE sources is almost irrelevant here since the y-ray flux is 
dominated by distant sources. We found that all models discussed in that sections reasonably 
reproduce the y-ray spectrum measured by Fermi-LAT and within errors agree with the model 
presented in M58I1 . We also found that the contribution of the extra-component electron to the 
bremsstrahlung and IC y-ray spectrum is subdominant and may be hard to detect. 

8. Conclusions 

The recent extension of the measurement of the CRE spectrum down to 7 GeV by the Fermi- 
LAT collaboration offers a strong evidence of the presence of a new spectral component. In fact, 
even allowing for a spectral break in the CRE source spectrum at few GeV, single component 
models fail to provide a satisfactory fit of those data. Such evidence adds to the one provided 
by PAMELA observations of a rising behavior of the positron fraction above 10 GeV. Assuming 
that the new component is symmetric in electron and positrons, we showed that a consistent 
interpretation of a wide sets of data is possible. For the first time, this also includes PAMELA 
positron fraction data even below 10 GeV as well as proton and antiproton data taken by the 
same experiment which are crucial to constrain the allowed propagation setup(s) and the solar 
modulation potential (during the same solar cycle phase Fermi-LAT is operating). We showed 
that only few among the propagation setups matching the B/C data can also consistently match 
CRE and positron fraction data. This is the case for Rraichnan-like and plain-diffusion setups 
while the commonly adopted Kolmogorov-type setup is disfavored. We obtained these results 
working in a simple force free, charge symmetric, solar modulation framework. This does not 
means that charge dependent effects, which are expected to be present at some level due to the 
complex structure of the solar magnetosphere, are absent but only that they seems not to be 
required for a consistent interpretation of Fermi-LAT and PAMELA data above few GeV. 

Concerning the origin of the e ± extra-component, we confirm that observed nearby pulsars 
are realistic source candidates. The expected anisotropy in the direction of the most prominent 
CRE candidate source, Monogem pulsar, is compatible with the present upper limit just released 
by Fermi-LAT collaboration and may be detectable in a few years. We also considered here the 
effect of known nearby SNRs on the CRE spectrum and the anisotropy, showing that they should 
not affect significantly those quantities. 

We conclude by observing that the AMS-02 space mission ll60ll to be installed on the Inter- 
national Space Station next year will be crucial to validate the models discussed in this paper and 
to look for small spectral anomalies which may indicate the presence of new physics. Beside 
providing and independent check of the positron anomaly detected by PAMELA and measure it 
up to 500 GeV or above, AMS-02 simultaneous measurements of CR nuclei, antiproton, electron 
and positron absolute spectra will allow to reduce the still large uncertainties on the propagation 
setup and the effects of solar modulation. Furthermore, its good energy resolution - which will 
be significantly better than Fermi-LAT - may allow to confirm the spectral features hints found 
in the Fermi-LAT CRE spectrum and to single out yet unidentified ones associated with nearby 
astrophysical sources or the dark matter. 
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